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Abstract 

We consider two chains of ultrasmall Josephson junctions, coupled capaci- 
tively with each other, and investigate the transport of particle-hole pairs and 
the quantum phase transitions at zero temperature. For appropriate param- 
eter ranges, the particle-hole pairs are found to play major roles in transport 
phenomena; condensation of such pairs leads to the superconducting state, 

displaying perfect drag of supercurrents along the two chains. 
PACS Numbers: 74.50. -Fr, 67.40.Db, 73.23.Hk 



Typeset using REVT^ 



Both the system of small metallic tunnel junctions and that of small Josephson junctions 
display the remarkable effects of Coulomb blockade, raising the possibility of single charge 
(electron or Cooper pair) tunneling [|I|,0. In particular, recent theoretical predictions 
and experimental demonstrations with two capacitively coupled one-dimensional (ID) 
arrays (i.e., chains) of submicron metallic tunnel junctions has revealed another fascinating 
role of the Coulomb interaction: The current fed through either of the chains induces a 
secondary current in the other chain. The primary and secondary currents are comparable 
in magnitude but opposite in direction. Such current drag at low bias voltage has the 
origin in the transport of electron-hole pairs, which are bound by the electrostatic energy of 
the coupling capacitance. Similar current drag has also been observed in two capacitively 
coupled two-dimensional (2D) electron gases 0, although the drag, attributed to a different 
mechanism [0], is much smaller in strength 0. On the other hand, in the Josephson-junction 
system, competition between the charging energy and the Josephson coupling energy is well 
known to bring quantum fluctuations, which affect in a crucial way the phase transition 
of the system |0-|T2|. Combined with such quantum fluctuation effects, the pair transport 
phenomena in coupled systems, which have not been studied, are expected to provide even 
richer physics. 

In this Letter we consider two chains of ultrasmall Josephson junctions, coupled capaci- 
tively with each other, and examine the transport behaviors and quantum phase transitions 
at zero temperature. For appropriate parameter ranges, the particle- hole pairs ( "excitons" ) , 
each consisting of an excess and a deficit in Cooper pairs across the two chains, are found to 
play major roles in transport phenomena. In particular, the condensation of such pairs leads 
to the superconducting state, which is characterized by the perfect drag of supercurrents 
along the two chains. 

We consider two chains of Josephson junctions (a = 1, 2), each of which is characterized 
by the Josephson coupling energy Ej [|T^ and the charging energies Elf'=e'^/2C\f' and 
E^'' =6^ / 2Ci'\ associated with the self-capacitance Clf' and the junction capacitance C{"'', 
respectively (see Fig. ||). The two chains are coupled with each other by capacitance C/, 
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with which the electrostatic energy Ei=e^ /2Ci is associated. It is assumed that there is no 
Cooper-pair tunnehng between the chains . 

It is convenient to write the partition function of the system in terms of the path integral 
representation 
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with the Euclidean action 
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The number n^°'\x) of excess Cooper pairs and the phase (f)^"'\x) of the superconducting 
order parameter on the grain at x in chain a are quantum mechanically conjugate vari- 
ables: [n^""^ (x) , (l)^''\x')]=iSxx'Sab- All the dynamics in the system occurs over the time scale 



-1- 



~-h/ y/8EjEs, where Eg = e^/2Cs with Cg determined by the largest capacitances (see 



below for precise determination). Accordingly, in Eq. the coupling constant has been 



defined to be K=y Ej/8Es, and all capacitances have been rescaled in units of Cg- Further, 
Vx and Vr denote the difference operators with respect to the position x and to the imagi- 
nary time r, respectively, and the (imaginary) time slice 6t has been chosen to be unity (in 
units of uj~^) [T^. The capacitance matrix in Eq. (|^) takes the form 



(3) 



where the submatrices C^"'\x,x') (a=l,2) are defined by the Fourier transforms 
C^''\q)=Cthc[''^ A{q) with A(g) = 2(l-cosg). 

Following the methods in Ref. |T^, one can transform the system in Eqs. (|1|) and (^ 
and, apart from the irrelevant spin- wave part, get the two coupled 2D systems of classical 
vortices: 
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where the 2D space-time vector notation r = {x, r) has been employed. The interaction 
?7af,(r) between vortices can be written in the block form 



C/(r) 



^ f/(2)(r)^ 
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where the subinteractions U^""^ and U^^^ are given by the Fourier transforms 
U(^){q) = [C^^\q)A{q) + D{q)A{uj)]/T{q) and f7(^)(q) = C/A(g)/T(q) with D{q) = 
C,C'(i)(g)+C(2)(g)] + C'(i)(g)C(2)(g) and T(q) = A\q) + /^(g)A2(^) + [2Ci + C(i)(g) + 
C(2)(g)]AHA(g). 

It is of interest to express the effective model in Eq. (0) in terms of the Hamiltonian 
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(a = 1,2,/) andt;(^)(r) = f W(r) - t;(2)(r). It should 
be noticed that for all Ci ^ 0, [/'^^-'(O) is always divergent and gives rise to the vortex 
number equality condition "^rV^^^r) = Z^r'^^^H'^)) equivalently, the vorticity neutrality 
condition J2r'^^^H''^) = Similar neutrality conditions Y.r'^^"'^ = (^^ = 1)2) should 



be satisfied unless Cq = [|T^. In general U^°'\r) for a = 1,2,/ is anisotropic in the x 



and the r directions, but the anisotropy is negligible at large length scales |jT2|; in the spirit 
of the renormalization group analysis, ^/'^"^(r) can therefore be considered isotropic in the 
space-time. In this way, the system has reduced to three subsystems of interacting vortices, 
as described by Eq. (|^): {v^^\v)} on (space-time) layer a (= 1,2) and additional vortices 
{v^^\r)}, which measure the relative displacements of the vortices on the two layers. Figure |^ 
shows that the configuration of the two vortices, f '•"'^^(r)=+l on layer 1 and f''^^(r')=+l on 
layer 2 (r 7^ r') gives a pair of vortex f '^■'^)(r)=+l and antivortex v^^\r')=—l for v^^\ For this 
reason, we call v^^^ "displacement vortices" . As we shall see in detail below, the displacement 



vortices describe effectively the correlations between the two chains in the system and play 
a major role in phase transitions and transport phenomena. 

The linear response ffabi^^) of the current in chain a to the voltage biasing chain b (see 
Fig. |I]) can be obtained via the analytic continuation 

o-a6(^) = — lim^ai,(g,iu;' uj+iO^), (7) 

where Qab is the Fourier transform of the imaginary-time Green's function 

Mx,r) = (T.[/('^)(x,r)/W(0,0)]), (8) 



with the time-ordered product and the current operators I^"'\x) = sin Va:0*-"''(a;). In 
the same manner as that leading to Eq. from Eq. (|I|), together with the path integral 
representation of the Green's function one can obtain the expression in terms of vortices: 

where (■ ■ ■)y stands for the average with respect to the effective vortex Hamiltonian in 
Eq. (D. 

Heretofore we have constructed a most general formulation for the system of two capaci- 
tively coupled Josephson-j unction chains: All the equilibrium and transport properties of the 
system can be understood through the 2D classical model in Eq. (|D or (||) and the response 
function given by Eqs. and (|^). We now consider a few simple cases and investigate the 
interesting phenomena associated with the particle-hole pairs. 

We first consider the case of two identical chains {Ci^^ = Ci^^ = Ci) without self- 
capacitance (Cq^"* = Cq^^ = 0)- III this case, with Cj = 0, each chain would be always 



insulating regardless of the ratio of Ei to Ej [O]. For large coupling capacitance Cj ^ Ci, 



the relevant capacitance and time scales are given by Cs = 2Ci and u;'^ = sJh^/4EjEi, 
respectively. At large (space-time) length scales, the short-distance anisotropy may be dis- 
regarded, which yields the vortex Hamiltonian 



r 



-^^eff E ^^^^r) log |r - r>(^)(r') (10) 



with the effective couphng constants k'^/J = ^Ej/16Ej and K^'^J = K^g = {Ci/Ci)K^^g . It 
should be noticed from the first term in Eq. ([T0|) that vortices on each layer does not even 
satisfy the vorticity neutrality condition and would form a plasma of free vortices. The state 
of the displacement vortices v^^\ on the other hand, depends strongly on K'^^g'- For k'^^^ 
below the Berezinskii-Kosterlitz-Thouless (BKT) transition value Kbkt^'^/'^ ||16[, they can 
exist in free vortices; for K^^g > Kbkt, only in bound dipoles. 

In general the system is insulating with free vortices (induced by quantum fluctua- 
tions) QJl^. Namely, the Coulomb blockade, disallowing single Cooper pair transport, 



tends to make each separate chain insulating. Nonetheless another transport mechanism 
is available for K^g > Kbkt'- To see this, we first note that the displacement vortices 
v'^^^ are nothing but topological singularities in the 2D (space-time) configuration of the 
phase differences (p'^^^ = (f)^^^ — 0*^^^ The absence of free vortices for K^g > Kbkt leads to 
the algebraic order of (p^^^ and accordingly the charge transport via the conjugate variables 
n^^^ = (n^^^—n^'^^)/2, each representing the excess Cooper pair number difference across the 
two chains or the number of excess- deficit pairs (i.e., particle- hole pairs). Such transport 
can be confirmed by examining the current responses given by Eq. (|^). Noting bound pairs 
of for Ki^ > Kbkt, one obtains at low frequencies from Eq. (0) 



TT 



Reanioo) ~ -Rea2i(^) ^ 2 - ^ 6{u), (11) 



where the renormalized coupling constant K^^ is given by 



1/4') = - y E |rP {v^'Kr)v^'K^))^^ . (12) 



r 



Thus the system exhibits superconductivity and carries currents along the two chains equally 
large in magnitude hut opposite in direction. This perfect drag of supercurrents reveals that 
the charges indeed transport in the form of particle-hole pairs, which are bound by the 
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electrostatic energy due to Cj. For K^^^ < Kbkt, on the other hand, the system displays 
insulating particle- hole IV characteristics, qualitatively the same as those in Refs. 
In a system with sufficiently small Ci {K^^g ^eff)^ therefore expect a BKT transition 
as Cj is increased, from the insulating {K^/J < Kbkt) to the superconducting state {K'^eff > 
Kbkt), where the charge transport occurs in the form of particle-hole pairs. The formation 
of bound dipoles of displacement vortices is an effective manifestation of the condensation 
of the particle-hole pairs, driving the superconductor-insulator transition in the system. 

One also observes similar behaviors in the system of two identical chains without junction 
capacitance {C[^^ = Ci'^ = and Cq'^ = C'q^ = Cq). Cq is assumed to be small so that 



both chains, with C/ = 0, would be in the insulating phase [0. We further assume that 



Cj ^ Co, and write the effective 2D vortex Hamiltonian in the form 

Hv ~ -Y.nKi^Y.^^''^^) log|r-r>('^)(r') 

- E ^^'^r) log 1 1" - r>(^)(r'), (13) 



where the effective coupling constants here are given by k'^Jij = K^J = yEj/SEj)'^ and 



k'^Jjj = ^Ej/16Ej. Unlike the previous case, vortices v^°'^ on each layer (a = 1, 2) as well as 
the displacement vortices v^^^ should satisfy the vorticity neutrality condition. Due to the 
condition i^^^^ <^ Kbkt, vortices on each layer separately would still form a neutral plasma 
of free vortices. Increasing Cj, however, we expect a BKT transition at k[^J = Kbkt, the 
superconducting state above which is again characterized by the condensation of particle- 
hole pairs and the current response function in Eq. (0). 

Finally, we consider the case Cq'^ ^ C/ ^ Cq^^ without junction capacitance {Ci'^ = 

(2) (2) 

Ci = 0), where the relevant capacitance scale is Cs=Cq . At large length scales, the 
vortex Hamiltonian is also given by Eq. (pTsD , with the effective coupling constants K^3f — 



Ej/8E^^\ K^^j^ ~ ^Cil&Q^ Kfjj, and k'^^J ~ K^J. We further assume that, 

with Ci = 0, chain 1 and chain 2 would be insulating and superconducting, respectively. 
This case is distinguished from the previous two in the fact that while chain 1 is insulating 
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[Re (Jii{uj) = 0], the response of chain 2 to the voltage biasing chain 1 is superconducting: 



R«^-H--;gy UM, (14) 



(2)- 



where Kj^'^ is given by Eq. (12) with the superscript / replaced by 2. This is not surprising 
since the particle-hole pair, with i?/ E'q , is no longer favorable. It should also be noticed 



that the current response in Eq. (|l^ is much smaller in magnitude than that to the voltage 
applied across chain 2 itself 

Rea22(^)^-(^ 2--f h(.;). (15) 

Therefore, although the particle-hole pairs still lead to the non-trivial current response of 
chain 2 to the voltage biasing chain 1, the quantitative effects are not so large as those in 
the previous cases. 

In conclusion, we have investigated the properties associated with particle-hole pairs in 
two capacitively coupled Josephson-j unction chains. In particular, the particle-hole pairs 
are found to play major roles in the transport phenomena and to drive the superconductor- 
insulator transition, with the superconducting state characterized by the perfect drag of 
supercurrents along the two chains. Although we have considered here the chains in the 



self-charging or nearest-neighbor charging limit, the argument in Ref. |12[ suggests that 
the quahtative results remain valid for realistic cases of generic capacitances. Such coupled 
chain systems can presumably be realized in experiment by current techniques, which have 
already made it possible to fabricate submicron metallic junction arrays with large inter- 
array capacitances as well as large arrays of ultrasmall Josephson junctions 0. We also 
point out that quasiparticles have been safely neglected in obtaining equilibrium properties 
at zero temperature. At finite temperatures or for large voltage bias, there may exist a 
significant amount of quasiparticles, which cause dissipation in the system [Q; still in the 
weak-dissipation limit, our results should not be affected qualitatively [|1^]. Further, it should 
be kept in mind that too large voltage or current biasing one chain can destroy the bound 
pairs of particles and holes, making the supercurrent in the other chain vanish. Finally, 
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note the difference between the straight couphng configuration 0, which has been studied 
here, and the slanted-couphng configuration [Q: In the case of tunnel junction arrays, the 
correlated sequential tunneling has been proposed as the relevant transport mechanism in 
the latter, whereas in the former the charge transport is attributed to the cotunneling of 
particle- hole pairs It would thus be of interest to compare the two configurations in 
coupled Josephson-junction systems. 

We are grateful to J. V. Jose for valuable discussions and sending us preprints. This 
work was supported in part by the BSRI Program, Ministry of Education of Korea and in 
part by the Korea Science and Engineering Foundation through the SRC Program. 
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FIG. 1. Schematic diagram of the system. 
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FIG. 2. Displacement vortex (a) The configuration of one vortex 'u(^)(r)=+l on layer 

1 and another v^'^\r')=+l on layer 2 (r 7^ r') corresponds to (b) a pair of displacement vortex 
■uU)(r)=+l and antivortex v^^\r')=—l. 
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